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Abstract 

The applied potential governs lithium-intercalation and electrode passivation reactions in lithium 
ion batteries, but are challenging to calibrate in condensed phase DFT calculations. In this work, 
the “anode potential” of charge-neutral lithium-intercalated graphite (LiCg) with oxidized edge 
planes is computed as a function of Li-content (uLi) at edge planes, using ab initio molecular dy¬ 
namics (AIMD), a previously introduced Li^ transfer free energy method, and the experimental 
Li^/Li(s) value as reference. The voltage assignments are corroborated using explicit electron 
transfer from fluoroethylene carbonate radical anion markers. PF^ is shown to decompose electro- 
chemically (i.e., not just thermally) at low potentials imposed by our voltage calibration technique. 
We demonstrate that excess electrons reside in localized states-in-the-gap in the organic carbonate 
liquid region, which is not semiconductor-like (band-state-like) as widely assumed in the literature. 


keywords: lithium ion batteries; voltage prediction; density functional theory; ab initio molecule 
dynamics; computational electrochemistry 
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I. INTRODUCTION 


The applied potential governs the thermodynamics and kinetics of lithium ion battery 
(LIB) interfacial processes. Li’*' insertion into graphite anodes to form LiCg is completed at 
0.1 V vs. Li’''/Li(s), or 2.9 V vs. standard hydrogen potential. Before this low voltage is 
reached, commonly used battery electrolytes containing ethylene carbonate (EC), cosolvents, 
and Li+/PFg salt already decompose at 0.7-0.8 V vs. Li+/Li(s). Fortunately, the growth of 
self-limiting films (called “solid-electrolyte interphase,” or SEI) formed via electron-injection- 
induced sacrificial electrolyte degradation passivates and stabilizes the anode.-i^ SEI con¬ 
siderations are also relevant for new anode materials like silicon. To intepret measurements 
and to help devise better artificial SEI/passivation layers,- there is a need to use electronic 
structure computational tools (e.g.. Density Functional Theory, DFT) to predict the voltage 
dependence of liquid-solid interfacial processes. Such a capability will have significant im¬ 
pact for studying not just LIB,-*^ but also lithium-air batteries,- water-splitting processes,-*^ 
and broad areas relevant to fuel cells, catalysis,-^— and electrodeposition.— In this work, 
we validate a recently devised potential calibration scheme,— apply it to interfaces between 
liquid EC and oxidized edge planes of LiCg, and explore the possibility of electrochemical 
decomposition of the counter-ion (PFg) used in commercial LIB electrolytes. In the elec¬ 
tronic supporting information (S.L), we document the signihcant error that can arise if the 
liquid electolyte is omitted in voltage estimates at interfaces. 

For non-redox active systems like non-Faradaic supercapacitors, the voltage difference be¬ 
tween two electrodes arises from their different surface charges mediated by electric double 
layers in the liquid regino.— In contrast, on complex LIB electrode surfaces, what a certain 
“applied voltage” means at the atomic level has not been sufficiently conceptualized, partly 
due to the difficulty in probing details at such lengthscales.— On the theory side, potential 
calibration has been challenging in periodic boundary conditions, condensed-phase DFT sim¬ 
ulations that depict liquid-solid interfaces.-'— DFT calculations are performed at constant 
number of electrons, not constant voltage. Each DFT simulation cell is associated with one 
electrode/Fermi level and is incompatible with a second reference electrode. Furthermore, 
most DFT electrochemistry calculations are conducted at T=0 K, which precludes explicit 
treatment of liquid solvents, dissolved salts, and, in the majority of cases, charged interfaces. 
Several interesting recent advances have coupled the Poisson-Boltzmann equation or related 
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approaches to DFT at the expense of introducing vacuum-solid or -liquid interfaces into the 
model.-'^ (The vacuum layer in effect serves as the reference electrode.) Purely condensed 
phase simulations with no vacuum region require other methods. 

We recently estimated what will be called the “anode potential” (V) of inert LiCg basal 
planes at hnite temperature.—^ The justihcation is briefly and heuristically described here. 
Our approach seeks to mimic LIB experimental processes, where Li-deintercalation occurs 
via transfer of Li"*" from LiCg, through the liquid-solid interface, to the liquid electrolyte, 
and ultimately into the Li metal counter electrode not explicitly depicted in the simulation. 
e~ flows in the same direction, but through the external circuit. Experimentally, it is known 
that these charge transfer processes occur at the onset potential of 0.1 V vs. Li’^/Li(s). 

We model this half-cell reaction at the onset of LiCg delithiation, 

(LiC6)n (Li(i_i/„)C6)„ -|-Li^(solv). (1) 

While our model does not include counter electrodes, and the excess e~ is left on the anode in 
the simulation (with hnite surface area correction, see below), Eq. [T]has effectively completed 
the e~ circuit. This is because, at equilibrium and in the absence of load-induced voltage 
drop, the Fermi level {E-p) of the Li metal “counter electrode” in our thought experiment 
must be lowered by 0.1 V to coincide with the Ep of LiCg. Under these conditions, the 
excess e~ on Lii_ 5 C 6 can start to how to Li(s). By reaching equilibrium (tuning the free 
energy change (AG^) of the Eq. [T]to zero) via varying the surface charge (a) on the electrode 
surface, we arrive at the experimentally known half-cell voltage for this reaction. This is the 
reference point that can be pegged to measurements. Away from this V=0.1 V vs. Li^/Li(s) 
hxed point, anode voltages are clearly related to the free energy of monovalent Li+ transfer 
between LiCg and EC liquid (AGt), AV = — AGt/|e|, provided that the interior Li atoms 
are frozen and not allowed to leak into the electrolyte (i.e., there is no redox reaction), which 
is the case in our simulations.— The a-V relationship associated with our frozen-Li basal 
plane LiCg model electrode is reminiscent of those in non-Faradaic supercapacitors. 

The ab initio molecular dynamics (AIMD, or DFT/MD) technique is used to calculate 
AGi.— AGt calculations deal with physical ions, not inhnitesimal/theoretical test charges, 
and therefore circumvent formal/unmeasurable concepts like the difference between “Volta” 
and “Galvani” dehnitions of the potential.-ii^ Further justihcation, including a thought 
experiment on explicitly including the interface between the liquid electrolyte and the Li 


3 


metal reference electrode, and comparison between our approach and related methods found 
in the aqueous computational electrochemistry literature, are given in Ref. |l5. 

The present work focuses on LiCe edge planes through which Li+ can intercalate and 
deintercalate. Edge planes are far more technologically relevant and complex than the proof- 
of-principle pristine basal plane considered previously.— But the same theoretical method 
can be used to examine the edge plane voltage. Indeed, Ref. (Fig. 2) implies that this 
Li+ transfer protocol can in priciple directly compare the voltage on any two electrodes. 
This is because a sufficiently thick liquid electrolyte region intervening between a basal and 
an edge plane can chemically and electrostatically screen them from each other. So the 
difference in equilibrium free energies of Li^ transfer from each electrodes to the electrolyte, 
divided by |e| should be proportional to their voltage difference. Since LiCg is an electronic 
conductor, at equilibrium e~ flows between surfaces to make the voltage the same on basal 
and edge planes. On pristine edge planes, the all-important AGi=0 point is determined not 
only by the electronic charge (a) compensated by mobile Li’*' in the liquid electrolyte, but 
also the fraction of Li occupying the edge sites (uli). At a hxed, applied potential, these 
two quantities should adjust themselves to minimize the free energy of the system. For 
the illustrative purpose of this work, we have hxed <7=0 and only varied uli, and V can be 
considered the instantaneous voltage before the edge Li’*’ content can change. 

Our predicted edge plane voltage as uli varies will be corroborated using explicit electron 
transfer from huoroethylene carbonate radical anion markers. Using AIMD simulations with 
appropriately calibrated potentials, PFg is shown to exhibit concerted e~ transfer and bond¬ 
breaking reactions at low voltages, suggesting that under such conditions electrochemical 
decomposition may need to be considered. This is signihcant because it is widely accepted 
that PFjT decomposes only thermally or due to reaction with trace water.— Finally, the 
simulation cells used in this work provide information about the electronic orbital alignment 
at explicit electrode/electrolyte interfaces. We show that excess electrons reside in localized 
states-in-the-gap in the organic carbonate liquid region, not in delocalized states at the 
conduction band minimum. The band structure is not semiconductor-like (band-state-like) 
as widely assumed in the literature.— 
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nu 

A(Li) 

^ ^max 

{dH{X)/dX)x AGt 

A 0.417 

0 

0.211 43.9 

-h5.45±0.11 

B 0.417 

0 

0.789 44.4 

-6.63±0.04 -1.14 

C 0.500 

0 

0.211 52.1 

-F5.86±0.10 

D 0.500 

0 

0.789 40.8 

-6.31T0.10 -0.78 

E 0.583 

0 

0.211 26.1 

-F6.59±0.11 

F 0.583 

0 

0.789 20.9 

-6.23±0.11 -0.36 

G 0.417 

1 

0.211 25.0 

-h5.80±0.05 

H 0.417 

1 

0.789 25.2 

-6.61±0.07 -0.94 


TABLE I: Details of AIMD trajectories for AGf calculations. A^(Li) is the number of mobile Li"^ 
in the liquid region, and A is the net charge of the Li^'*' ion frozen in the middle of the liquid, tmax 
is the total trajectory duration in picoseconds, and include the first 1 ps equilibration discarded 
when collecting statistics. The exceptions are C &: D, where tmax includes the discarded first 10 ps. 
Integrands and AGt are in eV; the latter is obtained by averaging the two integrands, and includes 
a —0.39 eV entropic correction and a —0.15 eV correction for using a 2-point treatment of Li'*’ 
solvation. 1^ To convert {—AGt) to V, add 0.1 V for Li'''/Li(s) and 0.1 V for the “half e~ rule” (see 
text). 

II. METHOD 


Our simplified electrode model consists of a LiCg strip with all C=0 termination.— >2^ 
Four neighboring C=0 bonds form a pocket where Li'*' can reside. The Li'*' surface density 
(nui) is unity if all such pockets are occupied. Figure [D depicts the periodically replicated 
29.74x14.97x15.06 interfacial simulation cell containing a Li3;Ci92048 anode slab and 
32 EC molecules. Interior atoms in the anode, with LiCe stochiometry, are frozen while C, 
O, and Li atoms at the edges are allowed to move, riu < 1 electrode conhgurations are 
obtained by randomly, and as uniformly as possible, removing Li atoms from the two LiCg 
edge planes. 

AIMD trajectories and a two-point thermodynamic integration (T.I.) formula with correc¬ 
tions/extrapolations are used to compute AGt (Table [11). These simulations apply the VASP 
code2^ with PAW pseudopotentials^^ and the DFT/PBE functional.— An energy cut-off of 
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400 eV, 10“® eV wavefunction convergence, and F-point sampling are enforced. Spot checks 
show that 1x2x2 fc-point sampling changes the integrands in AGt by less than ~0.05 eV, 
similar to basal plane cases.— A thermostat keeps the trajectories at an elevated T=450 K to 
improve statistics and prevent EC crystallization. Tritium masses are substituted for proton 
masses to permit 1 fs time steps. Compared to the previous work,^^ the predicted AGt is 
shifted by —0.17 V to correct the prior neglect of quantum nuclear effect inside bulk LiCe 
and the inadvertent use of T=450 K when adding translational and vibrational entropies to 
compare with experiments performed at T=300 K. 

The trajectories are initiated using Monte Carlo (MC) simulations in which anode atoms 
are frozen in DFT/PBE-optimized conhgurations. LIB salt concentration is typically 1.0 M, 
the static dielectric constant is large (Debye length ~3 A), and electrode surfaces should be 
screened from each other even in a small simulation cell. Electrical double layers should be 
well-equilibrated to the extent that the simple classical force helds used are accurate. More 
MC details are described in the S.I. 


III. RESULTS 


A. Controlling Potential at LiCg Edge Planes 


Instead of mapping the entire two-dimensional potential V{nu,o'), we focus on a = 0. 
In Fig. [1^, the linearly extrapolated V{a = 0,?7,Li) reaches the LiCe experimental plateau 
voltage of 0.1 V vs. Li'*"/Li(s) at nLi~0.69. If nLi>0.69, cr > 0 would be needed to raise 
V{a, flu) back to the green line and achieve the experimental potential associated with LiCe. 
This merely means that some of the edge Li must then be considered Li+ ions — not atoms 
— compensated with mobile PFg further away in the electrolyte.— 

Note that a “half-electron rule” vertical shift has been included to convert AGt (Table H]) 
to V(cT = 0,nLi) (Fig. [1^). T.I. calculations involve moving a Li+ from LiCg to the middle 
of the solvent region of a charge-neutral simulation cell, leaving an e~ behind. Along the 
T.I. path, an average of half an e~, or cr=—|e|/(4A), exists on the two electrode surfaces, 
where A is the lateral surface area. We assume this excess charge is uniformly distributed 
on the conducting electrode surfaces, in accordance with classical electrostatic predictions; 
see the S.I. of Ref. [iSl for analysis of instantaneous charge distributions. The effect of this 
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average a is estimated by finite difference and snbtracted from Fig. [T^, as follows. We inject 
one mobile Li+, compensating an excess e~ on the electrode, and recompute AGt (Tabled 
trajectories G-H). The simulation cells remain charge-neutral. AV is found to be —0.20 V 
(5cr/5V=17.8 /iC/(cm^G)) after adding the one e~. This is smaller in magnitude than that 
in basal plane simulation cells.To undo this surface charging effect, a -1-0.1 V correction 
is thus applied. The shift vanishes at large A, and represents an extrapolation to infinite 
system size. We stress that experimentalists can impose a potential without knowing details 
about the surfaces, but DFT calculations work differently; a and nu need to be adjusted to 
arrive at the desired voltage. 

In the absence of the liquid electrolyte, the potential at zero surface charge is directly 
related to the work function of the electrode in vacuum. In the S.I., we report work functions 
and show that the potential in vacuum, predicted as a function of nu, is significantly modified 
by the inclusion of the liquid electrolyte in the main text. This observation dovetails with 
predictions that work functions of metals can vary by ~1 V when their surfaces are covered 
with a monolayer of water— or organic molecules.— 

Although linear extrapolation is not expected to hold over the entire nu range, V(nLi, <J = 
0) appears to extrapolate to a value substantially below 0 V vs. Li+/Li(s) at nLi=l. Nega¬ 
tive potentials relative to Li~''/Li(s) are below the operating conditions of LIB anodes. This 
strongly suggests that a=0 and nLi=l yield an overpotential for electrolyte decomposition. 
In the literature, C=0 edge AIMD simulations have been reported at nLi=l and ct= 0, and 
EC molecules are found to decompose in picosecond via two different 2-e“ mechanisms, 
releasing CO and C 2 II 4 gases, respectively.— What are the potential dependences of these 
two competing processes? Recently, it has been predicted that the 2-e“ CO-releasing route 
has a far lower reaction barrier than C 2 H 4 generation in bulk liquid electrolyte regions.— It 
is therefore unlikely that C 2 H 4 should be a dominant 2 -e“ product unless there is an over¬ 


potential. Consistent with this deduction. Fig. [U indeed suggests the Ref. l22| system, where 


nu = 1, is at V < 0 V. If this were not the case, mostly CO products are expected. Note 
that two-electron reduction can yield C 2 II 4 gas if e~ are added sequentially, not simulta¬ 
neously, separated by milliseconds.— This timescale is far larger than our AIMD trajectory 
durations. 
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FIG. 1: (a) V at zero electronic charge as a function of edge plane Li content. The dashed line is 
a linear extrapolation, (b) A snapshot of the edge plane interfacial simulation cell. C, O, H, and 
Li are grey, red, white, and blue, respectively. 


B. Validating Predicted Potential: Electron Transfer 

The predicted potential should not depend on whether Li’*' or e~ moves across the inter¬ 
face. Next we demonstrate that V(cr = 0,?7,Li), calibrated using Li+ transfer above, is also 
consistent with e~ transfer. In the middle of the liquid region is placed a fluoroethylene 
carbonate (FEC“) radical anion (Fig. [2]), which is an effective electrolyte additive molecule 
for improving SEI on anode surfaces.— In charge-neutral FEC (Fig. [2]3), the carbonyl carbon 
(Cc) is coplanar with the three O atoms. In contrast, FEC“ is bent (Fig. [2]i), with Cc now 
sp^ hybridized. This leads to a large “reorganization energy” in the Marcus theory sense,— 
discussed below. We define the Cc out-of-plane displacement R as the scalar product be¬ 
tween (i) the normalized vector product connecting the three 0-atoms, and (ii) (Rc—Roave), 

8 








1 

w 




(e) 



FIG. 2: (a)-(c) i? as a function of time at various anode potentials V. Different line shapes denote 
trajectories with different initial configurations. (d)&:(e) Optimized FEC and FEC“ molecules. 
Out-of-plane coordinates R are illustrated. Purple circles represent F atoms. 

where “Gave” is the mean position of the three oxygen atoms. i?=0.0 A and 0.4 A in the 
optimized FEC and FEC“ geometries. Monitoring R grants easy access to FEC charge 
states. To initiate AIMD simnlations, we hrst condnct MC nsing classical force helds (see 
the S.I.), then rnn AIMD for 1.0 ps while freezing all FEC“ atoms, and hnally remove FEC 
constraints at “f=0” of AIMD trajectories. The time evolntion of R are shown in Fig. |2l 
The rednction potential of FEC in bnlk liqnid regions is predicted to be <F=0.58 eV when 
it is not coordinated to Li+ (S.I.). If V(ct = 0,nLi)<$, the excess e~ should stay on FEC“, 
ultimately leading to a second reduction of FEC“ and rapid decomposition.— >2^ In the 
opposite case, the excess e~ should be transferred to the electrode. The S.I. shows that this 
expectation is always satished with different initial conhgurations on less electrochemically 
active basal plane surfaces. At edge planes, this test is successful, but less than 100 % of 
the time. When V~0.18 V (nLi=0.667), lower than <F=0.58 V, FEC“ persists in all trials 
until it absorbs a second e~ and decomposes (Fig. [2^). This dovetails with our expectation. 
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When V=0.56 V, very close to <h, FEC“ is stable for hundreds of femtoseconds until it 
decomposes in two out of four trials; in the remaining two cases FEC“ loses its electron 
(Fig. [2b). Given our statistical uncertainties, this 50/50 split in the outcome is reasonable. 
When <F<V=0.98 V, FEC“ should give up its excess e~. This is observed 3 out of 4 times 
(Fig.[2t). 

On the whole, Fig. |2] demonstrates that V calibrated using Li’*' transfer also correctly 
governs e~ transfer. The one glaring “error” at V=0.98 V is apparently due to overly rapid 
C-0 bond breaking in FEC“ predicted using the PBE functional. This occurs within 100 fs 
(arrow in Fig. [21::; the predicted MP2 barriers for breaking this bond in EC“ and FEC“ are 
consistent with slower reaction rates.— Afterwards, the FEC“ ring cannot be reformed 
in AIMD timescale even if the anode potential favors it. It is also worth pointing out that 
DFT/PBE allows unphysical splitting of an excess e~ between the electrode and the redox 
center, artihcially accelerating e~ transfer rate,— and may make Fig. [2] unduly sensitive 
to initial conhgurations. Instantaneous fluctuations in the potential experienced by redox 
centers are physical and real, but e~ transfer occurs over a hnite timescale that partially 
averages out the fluctuations. Despite this, the overall trend of V and d> correspondence is 
correctly predicted: as V decreases, FEC“ retains its excess e~ more readily. 

C. Excess Electrons Form Localized States in the Gap 

Figure [3] depicts spatially-demarcated DFT Kohn Sham spin-orbital levels at f=0 ps in 
two trajectories taken from Fig. [2b & c. In panel (a) (V=0.56 V), the highest occupied 
molecular orbital (HOMO) is localized on FEC“ — reminiscent of a polaron in solid state 
physics —just below the Fermi level (Ep). FEC“ is found to decompose from this conhgura- 
tion. In panel (b), H=0.98 V, higher than <F=0.58 V; the localized orbital resides above Ep, 
and FEC“ loses its excess e~ rapidly in this trajectory. Fig. [3] illustrates that, for organic 
carbonate-based liquid electrolytes, the thermodynamic onset of electrolyte electrochemical 
reduction and SEI formation occurs when E-p coincides with the organic carbonate local¬ 
ized orbital found in the gap between liquid HOMO and liquid LUMO (lowest unoccupied 
molecular obrital or liquid conduction band minimum).— SEI formation does not begin, 
thermodynamically speaking, when the electrolyte LUMO coincides with the anode Fermi 
level, as has been widely assumed in the literature.— The LUMO lies above the organic 
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carbonate localized orbital and exhibits substantial statistical fluctuations; however it may 
influence e~ transfer kinetics if electron transfers from the anode, through a substantial liq¬ 
uid layer, to FEC in the bulk liquid electrolyte region, instead of towards FEC that diffuses 
near the anode surface. 

The above discussion adopts the solid state physics language often used in the battery 
community. It is important and of great interest to reconcile our study with molecular 
electrochemistry terminology^ less often featured in battery studies. According to Mar¬ 
cus theory,— e~ injection into FEC is accompanied with reorganization (free) energies (A). 
For FEC“, A contains a large intramolecular component, and is not solely due to “outer 
shell” solvation effects. The “polaronic shift” of the HOMO of FEC“, from above the liquid 
conduction band edge if FEC“ were flat, to within the liquid gap due to FEC“ geometry 
change and dielectric solvation, is a non-trivial manifestation of this A. Quantitatively, the 
vertical electron affinity should be at a value A above the molecular reduction potential 
Unfortuately, it is difficult to compute A in simulation cells with electrodes which are 
electron conductors. For example, the liquid electrolyte LUMO at frozen liquid geometry 
tend to reside above the Fermi level. As a result, injecting an e“ to the system to calculate 
the vertical electron affinity immediately populates the Fermi level of the electrode, not 
the the electrolyte, unless constrained DFT methods are used in the electrolyte region. In 
this sense, the systems considered in this work may differ from electrodes with a significant 
band gap, like Ti02.— Thus all AIMD simulations in this work report adiabatic free energy 
changes and redox potentials, not vertical excitations that include A. To estimate A, we 


have applied a localized basis set and dielectric continuum approach similar to Ref. 


29 


We 


find that the total A for FEC is 3.2 eV. This is somewhat larger than that of the struc¬ 
turally similar EC.— Our predicted A can be compared with future optical measurements 
in organic solvents analogous to those in aqueous media.— Incidentally, the wide separation 
between the localized excess e~ orbital in FEC“ and the liquid electrolyte HOMO minimizes 
hybridization between the localized state and solvent orbitals, which has been shown to be 
important for accurate DFT treatment of anions.— >2^ 
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FIG. 3: Instantaneous Kohn-Sham orbital decomposed on to atoms at their ^-coordinates at t=0 ps 
in Fig. c. 0 < x < 13.4 A denotes the metallic LiCe region; outside that range resides the 
liquid electrolyte which has a wide band gap. The Fermi level is at F1=0.0 eV. Panels a &: b depict 
the majority and minority spin channels of on trajectory each in Fig. [2^ and Fig. [2)3. The excess e~ 
“polaron” orbital of FEC resides in the majority spin channel, near E-p, in the electrolyte region. 

D. PFfT Electrochemical Reduction 

Finally, this potential-calibration technique is used to investigate possible PFg reductive 
decompostion on anode surfaces as V varies. A recent AIMD study has predicted rapid 
reductive decomposition of PF^ into LiF and PF^“ fragments, n <4, at liquid EC/graphite 
edge interfaces.— Seemingly consistent with this prediction, simple cluster-based calculations 
suggests that the LiPFg “molecule” has a reduction potential of 1.46 V (S.I.), far above LiCg 
potentials during battery charging. However, such a PFg reduction signature has not been 
observed in cyclic voltametry. Experimentally, the SEI formed in LiPFg-based electrolytes 
is known to contain LiF from PF^j" breakdown, but it has been widely accepted that PFg 
decomposes thermally or due to reaction with trace water over a period of hours,— >22 not 
electrochemically in seconds. 

To reconcile these observations, we note that, unlike FEC reduction, e~ transfer to PFg 
occurs in concert with P-F bond breaking. This is reminiscent of alkyl halide reduction,— 


12 

























which is clearly voltage dependent. Previons PFg modeling work— has not specified its 
anode potential. In Fig. lU we apply AIMD potential-of-mean-force (AkF(i?)) techniques to 
estimate the free energy barrier (AG*) of PF^ decomposition at two different potentials. 

The trajectories include a PFg pre-equilibrated at each graphite edge, charge-balanced 
by Li’*' in the electrolyte. The 3-atom reaction coordinate is i?'=|Rp — Rf| — iRu ~ Rf|; P, 
F, and Li atoms are chosen such that they are in position to react at t=0 ps. Increasing R' 
is correlated with F“ transfer from P to an edge Li. Umbrella sampling potentials of the 
form A{R' — are enforced, where 5=4 to 10 eV and Ro span the range between the 

reactant and the transition state. See the S.I. for details. 

FigureHK shows that, at V=0.56 V (?7,Li=0.583), AG* is at least 0.9 eV. As soon as one P-F 
bond breaks completely, e~ is transferred, other F“’s detach from the P-atom spontaneously 
(Fig. HJi), and these irreversible steps render a quasi-equilibrium sampling of AlU(i?) in the 
barrier top region impossible. In contrast, at V=—0.21 V (nLi=0.75), the barrier appears 
not much higher than 0.2 eV. Fig. 0] indicates that PF^j" electrochemical reduction may 
occur during the initial stage of SEI formation if the anode is at sufficiently low potentials, 
although this process faces competition from solvent reductive decomposition. 

Conclusions 

In conclusion, we have calibrated the anode potential (V((T = 0,nLi)) of lithium interca¬ 
lated graphite edge planes at zero surface electronic charge (cr=0) as a function of the edge 
Li content (n-pi) by computing the free energy of Li’*' transfer between electrode and liquid 
electrolyte. The estimated V is shown to be reasonable by correlating with observed electron 
transfer from reduced fiuoroethylene carbonate (FEC“) radical anions inserted into the liq¬ 
uid region. Electrochemical reduction of PFg at the pristine edge plane is shown to be viable 
at low potentials and to exhibit potential-dependent kinetics. This reduction pathway may 
need to be considered during SEI formation, in addition to thermal/impurity water-induced 
PFg decomposition routes widely accepted in the literature. In the future, optimization of 
the free energy with respect to all surface parameters {Nu and a) will be performed, and 
our method will be used to study the dynamics of Li+ insertion into passivated anodes as a 
function of the applied potential. 
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FIG. 4: (a) Potentials-of-mean-force (PMF) of PFg electrochemical decomposition, as functions 
R' at estimated 0.56 V (red, nLi=0.583) and —0.21 V (blue, nLi=o. 75 o)) respectively. The solid 
line portions depict actual data. The barriers must reside in the dashed line regions because PMF 
calculations there lead to spontaneous PFg decomposition in picosecond timescale, (b) Reduced 
Li^:F“:PF^ fragment from static, localized basis-set calculations. (c)&: (d): AIMD snapshots in 
PMF trajectories at F=0.56 V, ii'~—0.21 and R'~-|-0.5A, respectively. The rightmost PFg has 
decomposed in panel (d); recall that periodic boundary conditions are used. 
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